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Abstract. 

In a recent article |arXiv:1108.5127j Park has shown that the four-state predator- 
prey model studied earlier in J. Stat. Mech, L05001 (2011) belongs to Directed 
Percolation (DP) universality class. It was claimed that predator density is not a 
reasonable order parameter, as there are many absorbing states; a suitably chosen 
order parameter shows DP critical behaviour. In this article, we argue that the 
configuration that does not have any predator is the only dynamically accessible 
absorbing configuration, and the predator density too settles to DP critical exponents 
after a long transient. 



PACS numbers: 64.60.ah, 64.60.-i, 64.60.De, 89.75.-k 
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Systems having absorbing configurations may undergo a non-equilibrium phase 
transition [1] from an active to an absorbing state. The critical behavior of these 
absorbing state phase transitions (APTs) [2] depends on the the symmetry of the order 
parameter and presence of additional conservation laws. It has been conjectured [3] that 
in absence of any special symmetry the APT belongs to the directed percolation (DP) 
universality class as long as the system has a single absorbing state. 

Since the coarse grained microscopic theory of DP, which is a birth-death-diffusion 
process, is based on a single component Reggeon field theory [1], critical behavior in 
presence of additional field is expected to alter the critical behavior. The additional field 
may bring in multiple absorbing states and/or additional conservation laws. Presence 
of multiple absorbing states may [5] or may not [6] affect the universality. Coupling of 
order parameter to a conserved field too lead to DP [7] or non-DP [8] critical behavior. 
The models of directed percolation with more than one species [9], which brings in 
additional coarse grained fields, has also been studied [10]. The predator-prey cellular 
automaton models [H] in higher dimension too shows an APT to an absorbing (extinct) 
state which belongs to DP-class. The role of additional fields in these models are not 
quite well understood. 

Recently we studied a predator-prey model[l2] on a (1 + 1) -dimensional lattice, 
where each lattice site is either vacant, occupied by a predator A, a prey B or both (one 
A and one B). In these four state predator-prey (4SPP) model growth of preys and 
death of predators occurs independently, whereas death of a prey is always followed by 
instant birth of a predator. Based on the numerical simulations and estimated critical 
exponents, we have suggested the possibility of a new universality class. In particular, 
the decay of clusters at the critical point was found to be distinctly different from those 
of DP. However, in a recent article Park [13] has suggested a different scenario. It was 
claimed that the predator density ps can not be taken as a order parameter as there are 
infinitely many absorbing states. The transition is found to be in DP class, when order 
parameter is chosen suitably. In this article, we show that although there are many 
absorbing states, only one of them {p^ = 1,Pb = 0) is dynamically accessible. In fact 
the order parameter ps which was showing an apparently new critical behavior, slowly 
crosses over to DP. 

For completeness, first let us define the model. On a one dimensional periodic 
lattice, each site is either vacant, or occupied by a single particle A (prey), or occupied by 
a single particle B (predator) or by both particles (co-existing A and B); correspondingly 
at each site i we have Si = 0,1,2,3. Alternatively one may describe the 4-state 
predator-prey (4SPP) model considering two separate branches, one for A and the 
other for B particles, where particles living in one branch can not move to the other. 
We consider only the asymmetric case of the model which evolves according to the 
following dynamics, where X denotes both presence and absence of particles in respective 
branches. 
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Although, it was not exphcitly mentioned in our earher report [12], numerical 
simulations of the model was carried out in using exactly the same state variables 
hi = 4sj + Sj+i mentioned in [13]. Clearly these bond variables 6j = 0, 1 . . . 15 follow a 
dynamical rules, 

4 ^5; 6 ^7; 12 ^13; 14 ^15; 9 ^10; 13 ^14 (1) 
10^0;11^1;14^4;15^5; (2) 

The neighbors are updated along with hi as follows. For dynamics ([1]) fc^+i is increased by 
4 whereas for dynamics ([2]) fej+i and bi_i are decreased by 8 and 2 respectively. Clearly a 
bond i is active when bi is either 4 or 6, or it is greater than 8. It was argued by Park [13] 
that the density of active bonds p, the correct order parameter (as there are infinitely 
many absorbing states) of the system, vanishes at the critical point pc = 0.15381 when 
q = 0.02 and r = 0.9. For the same parameter values, we had estimated earlier that 
predator density ps vanishes at pf = 0.1484. This raises a question that possibly in 
the region pf < p < Pc, the system falls into an absorbing state which has isolated Bs. 
However a configuration with one or more isolated B can be absorbing only when there 
are no As. Thus, at pc we must have pA = 0. It was evident from Fig. 2(a) of Ref. [12] 
that pa> yp. This indicates that pb must vanish at the same value of p, possibly at 
Pc = 0.15381, as estimated by Park [13] . 




8 



Figure 1. (a) Plots of pBit)^^ vs i for p = 0.15375,0.15381,0.15385 with q = 0.02, 
r — 0.9 and L = 2^°. Here we use DP critical exponent a ~ 0.1594. (b) Both p{t) and 
PB{t) asymptotically decay as i-o is^^ _ 0.15381. 

To check this we redo the Monte-Carlo simulation for larger system size L = 2^° 
and measured psit) up to t = 3 x 10^ MCS for different values of p, while keeping 
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q = 0.02 and r = 0.9. From Fig. [T]^a) it is evident that pB{t)t°' shows a saturation at 
Pc = 0.15381 for a = acp = 0.1594. Figure [U^b) shows log-scale plot of p{t) and psit) at 
Pc', it is evident that after a long transient psit) decays as with a = 0.1594, similar 
to p{t). Such a change in a to a lower value was appearing as a saturation in ps, leading 
to a lower estimate of the critical point. 




Figure 2. (a) Log-scale plots of pB{t) vs t at criticality {p = 0.15381, q — 0.02 and 
r = 0.9)for L = 2^, 2^, 2^°, 2", 2^2 and 2". The initial slope a = 0.194 crosses over to 
aop — 0.1594 as system size is increased, (b) The same data are collapsed according 
to Eq. Q by using DP critical exponents (3/i'± = 0.252, and z = 1.580, 

With this correct estimation of critical point, we proceed to calculate other critical 
exponents taking ps as the order parameter. As for finite system of size L, starting 
from a high density of predators, Pg{t, L) decays as indicating a scaling form 

p,{t,L) = L-^/''-g{t/Ln, (3) 

at the critical point, where z is the dynamical critical exponent. Thus, Pg{t)Ll^/'^^ for 
different values of L are expected to collapse to a single function when plotted against 
t/L^. Figure |2]^a) shows decay of Pg{t) for different L starting from the configuration 
Pa = 1 = Pb', clearly small systems show an effective exponent a = 0.194 which crosses 
over to aDP as the system size is increased. True finite size effect of the critical point 
sets in at a reasonably large L. The data collapse according to Eq. ([3]) is observed in 
Fig. El^b) where we use the DP exponents ^ = 0.252 and z = 1.580. 

In summary, although there are many absorbing configurations in 4SPP model, the 
numerical simulations suggests that only one of them [pA = ^,Pb = 0) is dynamically 
accessible. The critical behavior of the absorbing transition can be well described by 
taking ps as an order parameter. The earlier estimated values of the critical exponents 
vary slowly with system size and settle to DP values. These studies truly emphasize 
the complications, difficulties and danger associated in numerical determination of 
universality class, when the possibility of a long transient is not ruled out. 
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